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Abstract 

This paper considers a multigateway multibeam satellite system with multiple feeds per beam. In 
these systems, each gateway serves a set of beams (cluster) so that the overall data traffic is generated 
at different geographical areas. Full frequency reuse among beams is considered so that interference 
mitigation techniques are mandatory. Precisely, this paper aims at designing the precoding scheme which, 
in contrast to single gateway schemes, entails two main challenges. First, the precoding matrix shall 
be separated into feed groups assigned to each gateway. Second, complete channel state information 
(CSI) is required at each gateway, leading to a large communication overhead. In order to solve these 
problems, a design based on a regularized singular value block decomposition of the channel matrix 
is presented so that both inter-cluster (i.e. beams of different clusters) and intra-cluster (i.e. beams 
of the same cluster) interference is minimized. In addition, different gateway cooperative schemes are 
analysed in order to keep the inter-gateway communication low. Furthermore, the impact of the feeder 
link interference (i.e. interference between different feeder links) is analysed and it is shown both 
numerically and analytically that the system performance is reduced severally whenever this interference 
occurs even though precoding reverts this additional interference. Finally, numerical simulations are 
shown considering the latest fixed broadband communication standard DVB-S2X so that the quantized 
feedback effect is evaluated. The proposed precoding technique results to achieve a performance close 
to the single gateway operation even when the cooperation among gateways is low. 
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Index Terms 


Multibeam satellite systems, multigateway transmission, MIMO precoding. 


I. Introduction 


Whenever multibeam satellite systems increase their overall throughput, their feeder link (i.e. 
the communication link between the gateway and the satellite) communication shall increase 
proportionally. Indeed, since the feeder link aggregates all satellite traffic, the resulting required 
bandwidth is 



( 1 ) 


where N is the number of feed signals, i?user-iink the per-beam bandwidth and i?feeder-iink the 
feeder link required bandwidth. Considering that the number of feed signals is proportional to 
the number of beams (K) and 


N > K, 


( 2 ) 


it is evident that if the user available bandwidth is increased, the feeder link resources must be 
increased accordingly. 

Augmenting the user available bandwidth can be performed if an aggressive frequency reuse 
among beams is deployed. In that case, interference mitigation techniques become mandatory as 
adjacent beams suffer from a large multiuser interference. Similar reasoning occurs in cellular 
terrestrial systems when spectrum sharing is considered [1]. 

Interference mitigation techniques are divided into two main approaches. One approach is 
that the receiver detects and, posteriorly, subtracts the multiuser interference. Examples of these 
advanced receivers can be found in [2, 3]. Another approach is to revert the interference effect at 
the transmit side, leading to high reduction of the receiver complexity since single user detection 
can be performed. Under this context, the transmitter must precode the feed signals in order to 
mitigate in advance the multiuser interference [4]. However, precoding is sensitive to channel 
state information (CSI) quality which is difficult to obtain due to the feedback and quantization 
errors. In any case, deciding where to implement the interference mitigation techniques exceeds 
the scope of this paper and in the sequel we will focus on precoding techniques. 

Apart from the already mentioned limitation of precoding, multibeam satellite systems have 
also several technical challenges [5]. As mentioned at the beginning, one of them is the feeder 
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link bandwidth requirements. In [6] the authors propose a feeder link system operating at the 
Q/V, where there is a larger available bandwidth than in the Ka band (eurrently used for this 
serviee). Although this frequeney band ean support the required feeder link bandwidth, rain 
attenuation severally deereases the performanee with respeet to the Ka band. Another option is 
to use multiple gateways and reuse among them the Ka band available bandwidth. As a result, 
the gateways must support the feeder links with very direetive antennas so that the bandwidth ean 
be fully reused among them. In any ease, whenever multiple gateways are employed, preeoding 
teehniques shall be reeonsidered [7, 8]. 

First, sinee the data traffie is independently generated at eaeh gateway, every gateway must 
preeode the signal in a deeentralized fashion and transmit through their eorresponding feeder 
link. In other words, the overall preeoding matrix is distributively eomputed at eaeh gateway so 
that eaeh gateway ean only use eertain feed elements. Seeond, sinee for eomputing the preeoding 
matrix eaeh gateway requires CSI from other gateways, eertain inter-gateway eommunieation is 
required, leading to large system resourees overhead. 

This paper deals with the problem of preeoding in multigateway multibeam satellite systems. In 
eontrast to previous works, this paper studies the ease where the eollaboration between different 
gateways is limited. In [9] a power optimization seheme is presented assuming linear mean square 
error (LMMSE) preeoding, full information sharing among gateways and single feed per beam 
arehiteeture. On the eontrary, our proposal eonsiders the multiple feed per beam strueture whieh 
is a more general payload arehiteeture and we study the ease where the information exehanged 
between gateways is redueed. 

Other works regarding multiple gateway preeoding are [8, 10]. In [8] a first approaeh to the 
problem is presented eonsidering that preeoding was able to bloek eompletely the interferenee 
between gateways whieh is not true in general due to the limited number of feed signals. 
Furthermore, in [10] a modifieation of the teehnique presented in [8] is deseribed eonsidering 
per beam power eonstraints. 

In addition to the aforementioned eontributions, in this paper some unexplored aspeets of 
multigateway multibeam preeoding are investigated. First, the tentative impaet of feeder link 
interferenee is eonsidered. Note that this effeet oeeurs whenever the payload feeder link reeeiver 
is not properly ealibrated and eertain pointing errors take plaee. Seeond, the limited eooperation 
among gateways is eonsidered and a matrix eompression teehnique is provided whieh shows 
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a good trade off between inereased overhead and overall throughput. Finally, the impaet of 
imperfeet CSI is evaluated eonsidering the eurrent standard digital video broadband (DVB-S2X). 

To summarize, the eontributions of this paper are: 

• A preeoding teehnique for multiple gateway arehiteeture where the transmit data and CSI 
is distributed among gateways. 

. The feeder link interferenee is studied and the mean square error (MSB) is eharaeterized 
under this effeet. 

• A proposal for redueing the eommunieation overhead between gateways is presented. 

. The effeet of limited feedbaek from the users and among gateways is evaluated assuming 
the latest standard for broadband satellite communieations. 

The rest of the paper is organized as follows. Seetion II deseribes the multigateway multibeam 
system arehiteeture so as its signal model. Seetion III provides a novel multigateway preeoding 
teehnique for multiple feed per beam payloads. In seetion IV is deseribed the limitations of 
multigateway multibeam preeoding teehniques. Seetion V shows the numerieal simulations and 
Seetion VI eoneludes. 

Notation: Throughout this paper, the following notations will be adopted. Boldfaee upperease 
letters denote matrices and boldface lowercase letters refer to column vectors. (.)^ and (.)^ 
denote Hermitian transpose and transpose matrices, respectively. Iat builds NxN identity matrix. 
Boldface 0 and 1 refer to an all-zero matrix and the vector of ones, respectively. Aij represents 
the (ith, jth) element of matrix A and {A)kxk denotes a submatrix of A of size K x K. diag 
represents a diagonal matrix. Tr denotes the trace operator. Finally, E{.} and ||.|| refer to the 
expected value operator and the Frobenius norm, respectively. 

II. System Description 

Consider a multibeam satellite communication system where a single geosynchronous satellite 
with multiple beams provides broadband services to a large set of users. The coverage area is 
divided into K beams and the users are assumed to be uniformly distributed within each beam. 

By employing a Time Division Multiplexed (TDM) scheme, at each time instant a total of 
K single antenna users, i.e. exactly one user per beam, is simultaneously served by multiple 
gateways. The total number of gateways is assumed to be G. In this context, each gateway 
serves a set of adjacent beams and the satellite is equipped with G feeder link receivers, each 
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one associated to a single gateway. Although in this first system model definition we do not 
include a tentative interference between feeder links, in subsequent sections we will investigate 
its possible impact. In any case, we will assume that the feeder link communication is ideal (i.e. 
noiseless communication) which is a realistic assumption considering the link budget of these 
systems. 

By assuming full frequency reuse pattern among beams, the received signals can be modelled 
as 

y = Hx + n, (3) 

where y G is a vector that contains the symbols received by each user, x G denotes 
the stack of the transmitted signals at all feeds and n G contains the stack of zero mean unit 
variance Additive White Gaussian Noise (AWGN) such that E{nn^} = Ik- Finally, H G 
is the overall user link channel matrix whose element (H)^ ^ presents the gain of the link between 
the i-th user (in the i-th beam) and the j-th satellite feed. Focusing on the channel model, the 
user link channel matrix H can be decomposed as follows: 

H = DW, (4) 


where: 

. D is assumed to he a K x K diagonal matrix which takes into account the atmospheric 
fading in the user link such that 

( ^ ^ 


D = diag 


( 5 ) 


V\/Ai '/Ar 
where denotes the rain attenuation affecting the k-th beam. 

• W is a if X matrix which models the feed radiation pattern, the path loss, the receive 
antenna gain and the noise power. The {k,n)-th entry of W is modelled as 

hVji gkn 


^kn 




( 6 ) 


where denotes the user receive antenna’s power gain, g^n is a complex value which 
models the radiation pattern from from the n-th to the k-th user, such that the respective 
feed transmit gain is 10 log^odfi'fenP) if expressed in dBi. Finally, dk is the distance between 
the k-th user and the satellite, A the carrier wavelength, ks the Boltzmann constant, Tr the 
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receiver noise temperature, BW the carrier bandwidth. Note that the elements of W are 
normalized by unit variance noise. 

The reader can refer to [11] for a more detailed description of the user link channel model. 
Assuming that the symbols are linearly processed before being transmitted, the transmit data 
can be written as 

X = Ts, (7) 


where s e is the transmit symbol vector such that the k-\h entry of s is the constellation 

symbol destined to the k-th user with E{ss^} = \k- T G denotes the block linear 

precoding matrix. Remarkably, T is computed at the gateways and it shall be transmitted through 
the feeder links. Thus, in the feeder link transmission the precoding matrix shall be separated in 
blocks and simultaneously transmitted by the gateways. This will be presented in the following 
sections. 

The typical performance metric involves the received signal to noise plus interference ratio 
(SINK) for k-th user, which is expressed as 

|(HT)fcfc|2 


SINRfr = 


Ef-I(HT),,|2 + 1 


k = 1, 


K. 


( 8 ) 


To generate a power flexibility which is essential for optimum resource allocation in multibeam 
system, we assume on-board travelling-wave tube amplifiers (TWTAs). In general, in TWTAs 
symbols are perfectly oversampled and pulse-shaped with a square-root raised cosine (SRRC) 
filter and a small roll-off factor. Therefore, these high power amplifiers provide a large bandwidth 
and high signal power level with small level of interference among feeds’ signals. In fact, the 
total power obtained by a set of TWTAs can be distributed among different feeds and the number 
of TWTAs is related to the number of feeds used in the antenna. Under this context, multiple 
feeds will distribute the available power to each beam. As a result, the precoding design is based 
on total power constraint. Mathematically, 

E{||x|| 2} = trace(TT^) < P, (9) 


where P denotes the total transmit power at the satellite. 

This paper is devoted to the study of multiple gateway precoding transmission so that in the 
following we will assume that beams are divided into clusters and each cluster is served by a 
certain gateway. There is the same number of gateways and clusters as Figure 1 depicts. As 
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we remarked previously, each feed signal can only be generated by a single gateway, otherwise 
signal overlapping might occur. Considering that there are G gateways, we will denote Ng the 
number of feeds associated to the ^'-th gateway for g = 1,..., G. Note that 

G 

JV = 5^JV,. (10) 

<?=! 

In addition, the number of served beams by the g-\h gateway is defined as Kg for g = 1, ... 
where 

G 

= (11) 

9=1 

In addition, we will consider that 

Kg<Ng g = l,...,G. (12) 

It is important to remark that the optimization of Kg and for o' = 1,... G is out of the scope 
of this paper and we will assume beforehand a certain feed and beam allocation per gateway. 
Consequently, each gateway must compute a precoding matrix G In other words, 

each gateway must compute the linear processing for serving Kg users with Ng feeds. Figure 1 
shows the overall system architecture. 

Prior to designing the precoding matrix, let us define the following division of the channel 
matrix 

H = (Hi,...,Hg), (13) 

where G is the channel sub-matrix containing the contribution of the feeds assigned to 

the g-th gateway. Note that we will assume that the allocated feed elements for each gateway are 
consecutive in the channel matrix. Moreover, each sub-matrix can be decomposed into clusters 
as 

= (14) 

where for c = 1,... ,G is the channel sub-matrix corresponding to the effect of 

the Ng feeds to the c-th cluster. An illustrative example of this sub-matrix decomposition is 
represented in figure 2. 

Next sections are devoted to the design of for = 1,..., G in order not only to mitigate 
the interference between clusters but also the interference created within each cluster. Posteriorly, 
limitations of this precoding technique are studied. 
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Fig. 1. Multiple gateway structure. The transmitted symbols are produced in geographically separated areas and they are 
transmitted separately through the feed signals. The satellite is equipped with an array fed reflector antenna where those N = 
^9 signals are transformed into K transmitted user signals. Interference is created not only by the signals within 
each cluster but also within the different clusters. The number of beams (users) per cluster is equal to Kg. 


III. Multigateway Multibeam Precoding 

As it is described in the previous section, the multiple gateway precoding does not only have 
to reject the intra-cluster interference but also the inter-cluster one. For the latter, it is proposed a 
block regularized singular value decomposition (SVD) precoding. Considering the ^'-th gateway 
channel submatrix, let us construct its regularized version 

Hf) = H,Hf + ^Ik e (15) 

which is used to construct the precoder instead of using Hg. The reason is that the null space 
has zero dimension and, as a result, it is impossible to block the inter-cluster interference. 
Nevertheless, as it will be shown in the simulation section, the out-of-cluster interference is 
mitigated severally if the regularized version is used instead. This technique was also applied in 
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Fig. 2. Representation an example of the channel matrix of a multiheam multigateway satellite system with three gateways. 
The grey suh-matrices represent the channel effect of the feed signals assigned to a certain gateway to the users assigned to the 
gateway. The white suh-matrices represent the impact of the feed signals to the non-intended clusters. The role of the precoding 
matrix is not only to revert the interference of the users within the same cluster (grey sub-matrix) but also to mitigate the 
interference generated to the rest (white sub-matrices). 


[12] for the multigroup multicast scenario. To understand how, let us decompose in the 
following sub-matrices 


H(«) = T:T(g,R),H 


(16) 


where G is the sub-channel matrix of the ^'-th cluster in the regularized matrix. 

In addition, the following matrix can be constructed with the out-of-cluster sub-matrices 




(17) 


Considering the singular value decomposition of this last matrix, 

= U.A^Vf, (18) 

where V^, can be written as 

= ( 19 ) 

where G are the left singular eigenvectors of associated to the non-zero 

singular eigenvalues, whereas V° G are the left singular eigenvectors associated to the 

zero singular eigenvalues. In this context, the gateway can use V° as a pre-processing in order 
to mitigate the interference created by its feeds to the users outside its corresponding cluster. 
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Now it is time to design the preeoding matrix in order to mitigate the intra-cluster interference. 
Considering that each gateway employs the aforementioned pre-processing matrix, the gateway 
observes a virtual channel matrix 

H^q = e (20) 

With this, the system designer can resort to the common designs such as zero forcing (ZF), 

= ( 21 ) 

Moreover, an additional gain can be obtained by means of employing the regularized zero forcing 

precoding or linear minimum mean square error (LMMSE) 

Wg^LMMSE = H?’" \ . (22) 

The final precoder is computed so that 

Tg = /5H^’^VjWg e (23) 

for both the ZF and MMSE case and where (3 is set so that 

Tr(TfT,)=| 9=1....,G, (24) 

where it has been assumed that the total power is equally divided for each gateway and the 

power constraints are fulfilled with equality. Finally, constructing the overall precoding matrix 
is done via 

T = block-diagonal{Tg}^,^^, (25) 

where the block-diagonal operator constructs a matrix formed by the inputs matrices in the main 
diagonal and zero otherwise. 

To sum up, the multigateway design consists of two stages. First, a pre-processing matrix that 
minimizes the inter-cluster interference is presented based on a regularized SVD decomposition. 
Second, for the resulting virtual cluster channel, two precoding techniques are considered: 
ZF and LMMSE. These latter techniques will mitigate the intra-cluster multiuser interference. 
Remarkably, the computation of the second precoding matrices can be done separately at each 
gateway where the first crucial pre-processing stage requires cooperation among gateways. 

In the next section precoding limitations are presented. Note that the multigateway precoding 
relays on the assumption that gateways know the overall channel matrix and; in addition, this 
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CSI has no errors. These two assumptions are generally unrealistic. The impact of limited CSI 
so as a possible reduction on the cooperation among gateways is studied in the next section. 
Finally, an analysis of the feeder link calibration errors is shown. 

IV. Multigateway Precoding Limitations 
A. Reduced cooperation among gateways 

The regularized block SVD multigateway precoding presented in the previous section assumes 
that each gateway knows H^. This assumption entails that each gateway is able to receive 
information by its feeder link from not only the users of its cluster but also from the rest. 
Unfortunately, this architecture will require complex data routing mechanisms at the payload 
which entails an unaffordable satellite hardware complexity. 

Assuming that the gateways are connected by a high speed fibre optic and assuming that 
satellite channel variations are minimal, full CSI sharing among gateways is implementable at 
expenses of a large communication overhead among gateways. In any case, the system designer 
might not decide to allocate such an overhead of resources for gateway cooperation and he/she 
can decide to keep them low. In that case, limited cooperation techniques are required. 
Concretely, the ^r-th gateway knowns 

g = l,...,G. (26) 

However, each gateway needs to know 

g = l,...,G. (27) 

As a result, each gateway should share 

{Hfigi,, 9 = (28) 

which is a total amount of Ng{K — Kg)G complex numbers to be transmitted through the 
connection between gateways. In order to reduce this communication overhead, the following 
approximation is presented. 

Considering the operation at the g-\h gateway, the CSI information to be transmitted to the 
i-th gateway is 

* ^ 9- 


g £KgXNg 


(29) 
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The best rank-one approximation of the aforementioned matrix is given by the right singular 
eigenveetor assoeiated to largest singular eigenvalue: cTjVg j G In this way, the ^'-th 

gateway ean eonstruet a matrix whieh approximates the inter-eluster interferenee as 

Gg = • • • ) . . . , ■ (30) 

Under this context, each gateway only needs to transmit Ng{G — 1) complex numbers instead 
of the Ng(K — Kg)G when full cooperation is considered. Matrix Gg E ^(G-l)xN, 

is an 

approximation of Similar to the previous method, the null space project of can be 
used as a pre-processing matrix following the scheme presented in the previous section. This 
development is not included in the paper for the sake of brevity. 

Apart from the presented reduction in terms of the singular value decomposition of Hf, the 
system designer can limit the cooperation overhead by means of reducing the communication 
between gateways. With this, only a subset of gateways interchange CSI data leading to a large 
reduction of the communication overhead. In the simulation section this is carefully evaluated. 

B. Feeder link interference 

So far we have considered that the gateways are able to transmit the feed signals within 
several feeder links (i.e. one feeder link per gateway) in an interference-free and noiseless 
channel. Indeed, the satellite is equipped with G feeder link receivers so that the signals are 
demultiplexed and routed through the array fed reflector. Unfortunately, although very directive 
antennas are used on ground for transmitting the feed signals, pointing errors might occur and; 
in addition, hardware might become uncalibrated under certain conditions. Consequently, an 
additional interference among beams might be created due to this. 

Mathematically, the channel matrix becomes 

(31) 

where H/- G models the interference generated by the different gateways and is the 

channel matrix described in the previous system model section. It is assumed that the interference 
equally impacts all the feed signals of each gateway (i.e the feeder link channel is not frequency 
selective). 
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In the following we provide a tentative deseription of Hj matrix. Assuming that is formed 
by bloek matriees, Yiy of dimension Ng x Ng whieh model the eaeh gateway interferenee, the 
feeder link interferenee impaet ean be deseribed as 


H 


/ 


9> 


hV = 


'-f ~ P' '"'^Ng * 7 ^ i) 


(32) 

(33) 


for z = 1,..., G and where Eat^ is a x Ng whose entries are equal to one. Moreover, p e [0,1] 
is a parameter that models the overall interferenee signals. The larger p the larger feeder link 
interferenee is eonsidered. 

Intuitively, whenever the systems presents larger interferenee prior the preeoding effeet, the 
lower aehievable rates will be obtained. This reasoning is mathematieally proved in the following 
by means of eonsidering an upper bound of the sum of MSB (SMSE). 

Let us eonsider that the inter-eluster interferenee is eompletely suppressed and eaeh gateway 
performs LMMSE preeoding. This seenario is an upper-bound of the overall multibeam multi¬ 
gateway system performance since in general the inter-cluster interference cannot be neglected. 
Under this context, the forward link system achievable rates under linear precoding considering 
an arbitrary g cluster can be described by the SMSE [4]. Mathematically, 


SMSE = trace (^MSe). (34) 

where MSE = diag^MSEi, ...,MSEj, ...jMSE^j and where MSEj refers to the MSE received 
by i-th user. Obviously, the lower SMSE, the larger achievable rates can be obtained. It can be 
shown that the SMSE can be written as 

SMSE = ^Trni^ + H„Hf^ V (35) 


where the g superscript has been omitted for the sake of clarity. With this last equation, the 
impact of the feeder link multigateway interference can be analytically studied. Consequently, 
it is possible to define 


SMSE 


no-interference — 



2l + H.Hf 



(36) 


K I [ G 

SMSEjnterference ~ Tr 11^ 




( 37 ) 
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where we did not include the superscript in H^- for the sake of clarity but it is important to 
remark that its dimensions are Ng x N. Prior to establishing the relation between the SMSE of 
these two cases, the following lemma is introduced. 

Lemma 1. Consider D = [di ,be any tall matrix and D,. = [di,...,dr], for all r = 
1,..., k — 1 it holds that 

cri(D^+i) > cri(D^) > (T2 (D^+i) > ... > ariBr+i) > cr^(D^) > (38) 

Proof. See Theorem 1 in [13]. □ 

With this result, the following theorem can be established. 

Theorem 1 For any semidefinite positive matrix H/, it results that 

SMSEjfitejfgj-ence ^ SMSEno-interference • (39) 

Proof See Appendix A. □ 

From this last result it is evident that whenever the multigateway feeder link structure is 
not perfectly calibrated, the system achievable rates decrease even though precoding is used. 
Thus, the multigateway feeder link hardware architecture shall be carefully designed in order to 
preclude the possible interference effect. In any case, as long as precoding is employed, the inter¬ 
feeder link interference can be reverted at the transmit side but; however, certain performance 
loss occurs with respect to the ideal interference-free feeder link system. 

C. CSI Feedback Errors 

Precoding severally relies on the CSI integrity. Preserving the quality of the estimation carried 
out by the receiver requires an ideal feedback mechanism which is impossible to implement in 
real systems. However, as it happens in terrestrial communications, broadband satellite standards 
are including CSI feedback mechanisms for supporting precoding techniques. This is the case 
of DVB-S2X [14] where for the first time certain feedback technique is offered to the system 
designer. In the following this mechanism is briefly described. 

The measurement and quantizing process of CSI feed back in each gateway is assumed to 
be continuous and to be reported on the return channels through a signalling table only when 
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significant changes are deteeted. The maximum delay required for estimation and delivery to 
the gateway via the interaetion ehannel shall be no more than 500 ms [14, Annex E.4], but this 
delay should be minimized to maximize eapaeity gain. 

Each user shall estimate and report the ehannel transfer funetion to the gateway as a set of 
eomplex-valued eoeffieients. These eoeflieients should be estimated by a set of 32 orthogonal 
Wash-Hadamard sequenees plus 4 padding symbols. In this context, every feed signal should 
ineorporate a different sequenee so that the reeeiver is able to estimate the ehannel effeet of 31 
interfering feed signals. With this, the overall ehannel matrix eannot be estimated but; however, 
the elosest 31 feed signals are the ones whose largest interference power levels so that the rest ean 
be ignored for preeoding purposes. Note that this also reduees the inter-gateway eommunieation 
overhead. 

Under this eontext, each user ean feed baek a maximum 7 digits (i.e. maximum 3 digits before 
deeimal point and 4 digits after deeimal point) for both amplitude and phase of eaeh ehannel 
element sueh that 

(H)^ . = ddd.dddd^'^"^^, (40) 

where it is evident that there are 7 bits for the magnitude so as 7 for the phase. Remarkably, this 
CSI report eonsiders the effeets of not only the user ehannel but also the tentative feeder-link 
imperfeetions. This is considered in the numerieal results section. 

V. Simulation Results 

A. Simulation Setup 

In order to illustrate the performanee of our proposal, this seetion presents the simulation 
results related to the eonsidered seenario in the previous seetions. The simulation setup is based 
on an array fed refieetor antenna whose ehannel gain matrix has been provided by ESA in the 
framework of a study on next generation multibeam satellite systems. The number of feeds and 
beams is assumed to be = 155 and K = 100, respeetively, whieh is eovering the whole 
Europe area such that N > K. 

We assume that the users are uniformly distributed over the coverage region and there is only 
one user per beam. Results have been averaged for a total of 500 user link ehannel realizations. 
Note that only atmospherie fading due to rain is considered in the user link ehannel and further 
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TABLE I 

USER LINK SIMULATION PARAMETERS 


Parameter 

Value 

Satellite height 

35786 km (geostationary) 

Satellite longitude, latitude 


Earth radius 

6378.137 Km 

Eeed radiation pattern 

Provided by if S'A 

Number of feeds N 

245 

User location distribution 

Uniformly distributed 

Carrier frequency 

20 GHz (Ka band) 

Atmospheric fading 

Just rain fading 

Roll-off factor 

0.25 

User antenna gain 

41.7 dBi 

clear sky gain 

17.68 dB/K 

Frequency 

20x10® 


refinements of the ehannel are not eonsidered. Furthermore, the randomness of the ehannel is due 
to the user positions whieh are assumed to be uniformly distributed within the beams. Reeall 
that, full frequeney reuse among beams have been eonsidered in this eontribution. Table III 
represents the detail of simulation parameters. 

We eompute the SINK for eaeh user, after interferenee mitigation and then its throughput 
(bit/s) is inferred aeeording to DVB-S2X standard for a paeket error rate (PER) of 10“® [14]. 
Table II provides a one-to-one relationship between the required reeeived SINR and the effieieney 
(bits/symbol) that are aehieved by the different adaptive modulation and eoding modes ineluded 
in the DVB-S2X standard. 

For a best praetiee and in order to elarify the performanee of proposed multiple gateway arehi- 
teeture, we eonsider the upper-bound of the system aehievable rates that result from eonsidering 
a single gateway system arehiteeture in [11]. 

B. Ideal Multigateway System Architecture 

This seetion presents the simulation results related to the seenarios deseribed in Seetion II 
and III. The user link in the downlink is assumed to operate at 20 GHz Ka-band. A total of 
G = 14 gateways are eonsidered so that eaeh gateway is serving a eluster of 7 or 8 beams. Fig. 
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TABLE II 

DVB-S2X MODCOD PARAMETERS 


ModCod 

Efficiency 

Required SINK [dB] 

mode 

Info bit / symbol 

(with approx, impl. losses) 

QPSK_2/9 

0.434 

-2.85 

QPSK_13/45 

0.567 

-2.03 

QPSK_9/20 

0.889 

0.22 

QPSK_11/20 

1.088 

1.45 

8APSK_5/9-L 

1.647 

4.73 

8APSK_26/45-L 

1.713 

5.13 

8PSK_23/36 

1.896 

6.12 

8PSK_25/36 

2.062 

7.02 

8PSK_I3/18 

2.145 

7.49 

16APSK_I/2-L 

1.972 

5.97 

16APSK_8/I5-L 

2.104 

6.55 

16APSK_5/9-L 

2.193 

6.84 

16APSK_26/45 

2.281 

7.51 

I6APSK_3/5 

2.370 

7.80 

16APSK_3/5-L 

2.370 

7.41 

16APSK_28/45 

2.458 

8.10 

16APSK_23/36 

2.524 

8.38 

16APSK_2/3-L 

2.635 

8.43 

16APSK_25/36 

2.745 

9.27 

16APSK_I3/18 

2.856 

9.71 

I6APSK_7/9 

3.077 

10.65 

16APSK_77/90 

3.386 

11.99 

32APSK_2/3-L 

3.289 

11.10 

32APSK_32/45 

3.510 

11.75 

32APSK_7/9 

3.841 

13.05 

64APSK_32/45-L 

4.206 

13.98 

64APSK_I1/15 

4.338 

14.81 

64APSK_7/9 

4.603 

15.47 

64APSK_4/5 

4.735 

15.87 

64APSK_5/6 

4.933 

16.55 

I28APSK_3/4 

5.163 

17.73 
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3 depicts the overall system architecture. It can be observed that the g-\h gateway is serving a 
set of 7 beams Kg = 7. Several collaborative schemes are presented in the following for the 
sake of completeness: 

• Scenario 1 ; The individual cluster multibeam processing without gateway coordination so 
that each gateway processes its beams independently. In this context, the ZF precoding in 
m-th gateway can be expressed as 

T, = AcMHf (H^Hf )-\ (41) 

where (5icm is set to preserve the gateway power constraint assumed to be ^ for all of them. 
This is referred to Individual Cluster Multibeam processing (ICM). 

. Scenario 2 : 4 gateways interchange their CSI so that each gateway only has access to 3 
interfering matrices. This is referred to 4 Gateways Collaboration (4GC). See Fig. 4. 

. Scenario 3 ; 7 gateways interchange their CSI so that each gateway only has access to 3 
interfering matrices. This is referred to 7 Gateways Collaboration (7GC). 

• Scenario 4 : Gateway g (respectively for all the gateways) collaborates with all gateways. 
This is referred to Gateway Collaboration Multibeam processing (GCM). 

. Scenario 5 ; Single gateway scenario where all data and CSI is located at the same 
transmitter. This scenario is the Reference scenario (Ref). 

. Scenario 6 : Gateway g (respectively for all the gateways) collaborates with all the gateways 
that serve clusters that are directly adjacent to cluster g by means of transmitting the rank 
one approximation of their channels as we discussed in section IV A. This is referred to 
Limited Multi-gateway Collaboration processing (LMC). 

Table III describes the communication overhead associated to each cooperative scheme. As 
we anticipated in the previous section, the LMC scheme offers a large reduction of the commu¬ 
nication overhead. 

Figure 5 and 6 present the spectral efficiency when either MMSE or ZF is used as precoding 
matrix for mitigating the intra-cluster interference respectively. For both cases it is shown that the 
ICM scheme has the lowest achievable rate due to lack of interference mitigation among clusters. 
4GC and 7GC achieve reasonable performance and the proposed multiple gateway scheme with 
CSI sharing among adjacent clusters (i.e. GCM) shows better performance. Consequently, the 
higher the coordination among gateways is, the higher the achievable rates are . Note that MMSE 
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TABLE III 

Cooperation Overhead Comparison 


Cooperation Scheme 

Total Number of Complex Numbers to be Shared 

4GC 

57288 

7GC 

100254 

GMC 

200508 

LMC 

143 



Longitude [Deg] 


Fig. 3. A set of 14 cluster composed by 7 beams is depicted in the figure. This will be the reference architecture assuming that 
all beams operate in the same frequency band. Remarkably, the precoder not only has to mitigate the inter-cluster interference 
but the intra-cluster one. 


delivers higher speetral effieieneies than ZF. 

Note that the proposed LMC offers a good trade-off between gateway eooperation overhead 
and overall system performanee. 

C. Non-Ideal Feeder Link System Architecture 

In this seetion we analyse the seenario presented in the seetion IV. We eonsider p = 1, whieh is 
a worst-ease seenario. Again, simulation results use the average total throughput as performanee 
measurement. Figure 7 eompares the results related to all seenarios deseribed above eonsidering 
that for eaeh gateway reeeives interferenee from 1 to 14 gateways. The transmit power is set to 
F = 30 dBW and MMSE preeoding is used for mitigating the intra-eluster interferenee. 

From Figure 7 it is evident the dramatieal effeet of feeder links mismatehes in multi-gateway 
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One cluster 


o 

J? 



Longitude[Deg] 


Fig. 4. Cooperation architecture between 4 adjacent clusters: it is assumed that each gateway can cooperate with only 3 gateways 
whose beams are adjacent to them. The cooperative clusters are depicted with the same colour. 



Fig. 5. Spectral efficiency considering multigateway block regularized precoding and different collaborative architectures. The 
intra-cluster interference is mitigated via MMSE precoding 
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Fig. 6. Spectral efficiency considering multigateway block regularized precoding and different collaborative architectures. The 
intra-cluster interference is mitigated via ZF precoding 


multibeam architecture even though precoding is performed. This effect appears for any co¬ 
operative architecture. It is important to remark that even if only one interfering gateway is 
considered, the average spectral efficiency decreases up to the 54% with respect to the ideal 
feeder-link scenario. 

D. Limited Feedback 

Let us consider the case where the CSI is not perfect but it is obtained via DVB-S2X feedback 
mechanisms as we described in section IV B. 

Figure 8 depicts the performance degradation whenever quantized feedback is used instead of 
the perfect one. Both for the LMC and GMC case, the spectral efficiency decrease is notorious. 

VI. Conclusion 

This paper presents a precoding design for multigateway multibeam satellite systems. These 
architectures suffer from a large intra and inter cluster interference. The proposed technique 
mitigate both interference types, leading to large spectral efficiencies. Although the presented 
method assumes that gateways share the CSI (i.e. all channel matrix), limited cooperative schemes 
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Fig. 7. Spectral efficiency versus number of interfering gateways. The interfering parameter is set to 1 (p = 1). The feeder 
links uncalibrations decrease the system capacity severally. 



Fig. 8. Spectral efficiency versus transmit power considering GMC and LMC processing with CSI errors. The performance 
decrease is high whenever the available CSI is not perfect. 
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are introduced. In addition, whenever certain feeder link imperfections are considered, it is shown 
that the SMSE and thus the achievable rates decrease. The proposed techniques were validated 
through the DVB-S2X standard which recently proposed certain CSI feedback mechanisms. 


Appendix A 
Proof of Theorem 1 

The SMSEs can be rewritten so that 

i<r 

1 


SMSE 


interference 


K 


1 


K 

SMSEj^o-interference 


E 




(42) 


(43) 


Under this context and manipulating the previous equations, the theorem 1 holds as long as 


A.(H,HfH„H") < Ai(H«H„), 
for z = 1,..., iT. Bearing in mind that, 

Az(H„H^HfH{) = 

Az(H„Hf) = ), 


(44) 

(45) 

(46) 


and considering that has the following SVD decomposition = U/SjVj?, we have that 


H;Hf = U/S;Sfvf. 


(47) 


Writing = S/ it is easy to observe 

'^NgXNg 0NX{N-Ng) 

0{N-Ng)xNg 0(N-Ng)x{N-Ng) 

where (S/S^) has only Ng non-zero singular values (i.e. z), as Sg has rank equal to N^- 
Since U is a unitary matrix, for any matrix (same as in this study) it holds that 

a,(UHfH„) =ai(HfH„). (49) 

Then, the right hand side of (44) can be worked out as 


S,= 


(48) 


a, (U^S^UfH^H.) = a, (S^UfHfH,,,) . 


( 50 ) 
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By the following definition Gf = U^H^Hu and eonsidering (49), proving (44) is equivalent 
to eheeking 

(51) 

for i = 1, Ng. Now, remind that G/ is of size N x N as follows 

^ (g:) 

where both sub-matriees of G/- are of size Ng x N. Then, we have that 

S/G/ = (53) 

It is elear that 



where G^ is matrix of size N x N whose singular values are 


cr(Gf) - (cri(G/),(T2(G /),... ,aNgiGf),ON-Ng) , 


(54) 


(55) 


where is a veetor of dimension N — Ng whose entries are equal to zero. By the interlaeing 

property in Lemma 1 tells that 

a^iGf) > a.iSfGf). (56) 
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